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Abstract 

We address the issue of stock market fluctuations within Langevin Dynamics (LD) and 
the thermodynamics definitions of multifractality in order to study its second-order char- 
acterization given by the analogous specific heat C q , where q is an analogous temperature 
relating the moments of the generating partition function for the financial data signals. 
Due to non-linear and additive noise terms within the LD, we found that C q can display 
a shoulder to the right of its main peak as also found in the S&P500 historical data 
which may resemble a classical phase transition at a critical point. 
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1 Introduction 



Simulated descriptions of at least some universal features of financial prices has recently 
been produced by multifractal techniques U and also linear O, |I| and non-linear |4j] 
Langevin Dynamics (LD). On one hand, multifractals reasonably describe patterns that 
resemble sudden (both large and small) market price oscillations, ensuing bursts of long- 
range volatility as in the case of stochastic multi-agent models ||. On the other hand, 
non-linear LD generates fat tail distributions of price differences as found in the 
analysis of financial time series || |7j . 

There is a vast literature on the characterization of financial data based on tur- 
bulence, multifractality, autocorrelation functions and power spectra studies (see, e.g., 
||, [|, |nj [11], |~2|). However the second-order derivative relating the increment of price 
values of assets Ax(t, 1) = x(t + 1) — x(t) and given by an analogous specific heat C q 
similarly to the multifractal studies of diffusion-limited aggregates (DLA) |L3| , has only 
recently been analysed by the author in the case of temporal fluctuations of the S&P500 
stock index ||14|| . The analogous specific heat is given by the second order derivative 
of the analogous free energy r relating the multifractal dimension D q (c.f., Eq.(|]) and 
(0) below). The g-moments of the generating partition function for the financial data 
signals Ax(t, 1) corresponds to an analogous temperature. 

As an illustrative example, the curve (a) in the lower Fig.l shows the presence of 
an anomalous shoulder in C q when considering 3287 (S&P500) data points including 
the largest burst that corresponds to the so-called Black Monday crash measured in 
October of 1987. For comparison, in the plotted curves (b,c,d) of Fig.l we remove the 
crash from the historical data and make the measurement of C q anew. Curve (b) is the 
result of just removing the single data point for the Black Monday (i.e., by considering 
3286 points), whereas the curve (c) is the result of removing the Black Monday plus 
all successive data {i.e., considering 1970 initial points). On the other hand, curve (d) 
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is the result of taking 499 initial data points only. 

Clearly, as one removes or approaches the Black Monday data the shoulder in C q 
of curve (a) does not disappear as can be seen in curves (b,c). A similar behaviour is 
observed analysing the data from few days prior to the S&P500 crash. As a difference 
C q only displays a typical, single peaked form without anomalies for the data prior to, 
but far away from, the Black Monday (as in curve d). 

In this work we make an attempt to characterize the onset of market crashes by 
analysing the second-order derivative relating the increment of prices of assets Ax(t, 1) 
and given by C q . We shall consider the LD of financial systems || [| to show that the 
presence of an anomalous shoulder in the specific heats as those plotted in Fig.l can 
also be a consequence of having non- linear stochastic terms within LD. 

2 Linear and Non-linear Langevin Dynamics 

We consider the differential equation with an additive noise term of the form g t times 
Vt- 

= h t (x(t)) + g t (x(t)) Vt , (1) 

where as usual rjt is assumed to be Gaussian-white whose average and variance are given 
by r}t = and rfirg = 2DS(t — t'), respectively. 

The dynamics of simplest financial systems with random noise has been proposed in 



I5| , |16fl by setting 

h t (x(t)) = -X t x(t) ; g t (x(t)) = l , (2) 

with X t a multiplicative noise that sets the time scale for deviations from equilibrium. In 
the context of economic models, x(t) may represent price values, A t the rate of change 
of the price, and r)t some external noise of various sources. In the following we shall 



refer Eq.(|2|) as the linear LD. We add that the scaling behaviour and the 1/f noise of 
trajectories and velocities with time within linear LD with fluctuating random forces 



has been reported in our previous work [ |17|j . 

Another simple financial agent model has been introduced in [f| in analogy with the 
mean field approach to the Ising model for a magnetic system characterized by 

h t (x(t)) = Jx{t) + bx 2 (t) -cx 3 (t) ; g t (x(t)) = e + x(t) , (3) 

with J, b and c constants associated with higher-order agent interactions leading to power 
laws tails. The parameter e allows this agent model to take up the simple Langevin form 
when x(t) — > 0. The dynamics of time series from stochastic processes governed by a 



similar non-linear LD has also been reported in |18 . 



3 Multifractality and Analogous Specific Heat 

In accordance with the standard of economic notation only relative changes Ax(t, 1) are 
usually relevant. From these fluctuations, we follow |T^] and construct the generating 
partition function Z in terms of the normalized measures fit > 0, and its moments q, by 
the scaling 

N 

Z( qi N) = Y.tf~N~^) , (4) 
t=i 

where 

_ |A.x(t,l)| 

^ = Ef =1 |AaKM)| • (5) 
In the above we have divided the ID system of length L into N lines of length I (i.e., 



N ~ L/£), and have associated this iV with the measured discrete x(t) sequences |L9 
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In order to get the thermodynamics interpretation of multifractality (see, e.g. |fL3[|), 
we then consider standard definitions 



T{q) = [q-1]D 9 , (6) 

where r represents an analogous free energy and D q the multifractal dimension. Of 
particular interest here is to consider the analogous specific heat given by the second 
order derivative of r, namely 

-C 9 = ^^«r(g + l)-2r(?)+r(?-l) ■ (7) 

4 Discussion 



In the upper curves of Fig. 2 we show some typical time evolution of Ax(t, 1) governed 
by the linear LD with h t and gt given in Eq.(^) for different values of rjt = (1 — 2r n >)v 
and A* = (1 — 2r n )u with v = 0.1, u = —0.005, —0.02 and —0.1, and using the random 
numbers < r n>n > < 1. In the lower curves of Fig.2 we display the non-linear LD 
fluctuations of Ax(t, 1) with h t and g t given in Eq.(||) with J = —1, b = 0.01, c = 0.001, 
e = 0.005 and for different values of rjt — (1 — 2r„»)n such that n = 1.2, 1.5 and 1.8 and 
using random numbers < r n » < 1. 

From the curves in Fig.2 it can be seen that Ax(t, 1) has relatively small finite 
values, and exhibits intermittent bursts due to the weak noise terms. The presence of 
peaked narrow fluctuations can increase with both the multiplicative and additive noise 
strength. However, wider volatility clusters are formed over large periods of time when 
decreasing, in the linear LD, the amplitude u of the multiplicative At noise as compared 
to the case of increasing n of the additive noise rj t within the non-linear LD. It is this 
difference of behaviour that will lead C q to display an anomalous shoulder to the right 
of its main peak to be shown below. 
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Let us first see if the generalized dimensions D q for both linear and non-linear LD, 
which are extracted by considering the definition in Eq.(^|) in conjunction with the scaling 
of Eq.([|), are sufficiently smooth and multifractal for C q of Eq.(^) to be meaningful 
|20| . |2~T| . The discrete values of q used range from —5 to +5 at increments of 0.033. 

In Fig. 3 we display the plots of D q for the different values of u and n as used in Fig.2 
together with the plots of the limit 1 — for the linear and non-linear LD (upper 

and lower curves, respectively). If q < — 1, we observe that both results for D q follow a 
typical convergent behaviour according to multifractal physics. In the range u > —0.02 
and n < 1.5, D q is fully multifractal-like. 

For lower values of u < —0.02 and greater values of n > 1.5, we find a non- 
monotonous decreasing behaviour of D q , in correspondence with the double-peaked 
form of the respective C q displayed in Fig.4 which relates the presence of the inter- 
mittent bursts in Ax(t, 1) shown in Fig.2. Such a non-monotonous behaviour of D q is 
more evident in the case of non-linear LD. From the D q data for the linear LD case, 
the multifractality strength of the Ax(t,T) = x(t + T) — x{t) sequences at different 
integer time lags T > 1, i.e, 1 — D^^, does not seem to follow a power-law scaling for 
1 < T < 10 within all u considered. As a consequence of the complicated behaviour 
of the non-linear form of Eq.(Q) coupled to i] t , which influences the (thinner) volatility 
clusters shown in Fig.2 (lower curves), it is not strightforward that 1 — i n this 

case follows a power scaling as in PH . 

Let us see next how this type of D q behaviour given in Fig. 3 influences the analogous 
specific heat of Eq. (0) and how it characterizes the temporal large and small intermittent 
bursts representing variations in financial prices. In Fig.4 we plot C q for the LD param- 
eter set used in Figs. 2 and 3. For the linear LD case, we find that the main peak of our 
numerical C q curves for u > —0.02 resembles a classical first-order phase transition at 
a critical point q = —1. However, the sharp peak turns slightly asymmetric around the 
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critical point for u = —0.1. Surprisingly, this asymmetry becomes evident within the 
non-linear LD. The analogous specific heat in this case displays a higher shoulder to the 
right of the main peak as a function of n > 1.5 which is in qualitative agreement with 
the C q curve for the S&P500 index data shown in Fig.l (lower curves). The presence 
of the shoulder in C q is then associated to the non-monotonous behaviour of D q as in 
Fig.3. 

The same C q shape at a phase transition as in Curve (d) of Fig.l and Fig. 4, has 
been reported in the case of DLA |i~3| . A non monotonic specific heat shape with a 
superimposed second peak (resembling the shoulder of the type of Curve (a) in Fig.l 
and Fig. 4 lower curves) has been found in the Hubbard model on a cluster of magnetic 



sites p2| , |23|j and, most recently, in a uniform spin model on a fractal | 23fl . 

When decreasing n < 1.2 in the non-linear LD, the sudden intermittent bursts in 
Ax(t, 1) as well as the shoulder in C q tend to vanish. It is this latter feature that make us 
argue that the stochastic price behaviour of financial assets may be characterized by an 
analogous C q which resembles the phase transition features as measured in multifractal 
physics. 

5 Conclusion 

We have made a novel attempt to characterize the onset of market crashes via an analo- 
gous specific heat using the thermodynamics definitions of multifractality. We presented 
an study of generalized dimensions D q of Eq.@, and the analogous specific heat C q of 
Eq.([?D, characterizing the absolute moments of the increments of two LD, a linear one 
(c.f., Eq.@) and a non-linear one (c.f., Eq.(|3])), within the context of a multifractal 
approach to financial time series analysis 



Tunning parameters in these dynamics introduces wide clusters and bursts of volatil- 
ity in the processes. This behaviour was shown to mimics characteristics of S&P500 
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data which were measured according to an analogous phase transition in the behaviour 
oiC q . 

In this work, the second-order derivative relating the increment of prices of assets 
Ax(t, 1), and given by C q , follows similar Legendre transforms for the thermodynamics 
analysis of (DLA) ||13|| . Our findings may relate the results of the kind of those described 
25| where a phase transition-like appears in the presence of two distinct phases when 



in 



analysing the behaviour of D q . In the latter, these two phases are built out of rare events 
on the one hand and some "background" (e.g., bursts of volatility in the processes) on 
the other hand. We have found that a shoulder appears in the analogous C q due to the 
presence of wide volatility clusters as found in the S&P500 historical data, or in the 
non-linear LD. 

The message here is that due to non-linear and additive noise terms within the LD, 
C q can display a shoulder to the right of its main peak data as also found in the S&P500 
historical data, which may resemble a classical phase transition at a critical point. 
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ure captions 



Fig.l: Time evolution Ax(t, 1) of the S&P500 stock index for the period 1980-1992 
as a function of trading time lags (in a.u.) and corresponding analogous specific 
heat C q obtained using Eq.(^). The lower curve (a) shows the presence of an 
anomalous shoulder in C q when considering 3287 (S&P500) data points including 
the largest burst that corresponds to the crash measured in October of 1987. Curve 

(b) is the result of removing this Black Monday crash data point, whereas curve 

(c) is the result of removing all successive data. Curve (d) is the result of taking 
499 initial data points only. 

Fig. 2: Time evolution of Ax(t, 1) governed by the LD of Eq. (Jl|) using the linear 
terms in Eq.(|2|) (upper curves) and non-linear terms of Eq.(|3]) (lower curves). 

Fig. 3: Generalized dimensions D q and multifractality strength 1 — -D^oo for differ- 
ent time langs T for the linear LD using Eq.(|2]) (upper curves) and the non-linear 
LD using Eq.(|3]) (lower curves). The 1 — D^^, curves have been evaluated by 
considering Ax(t, T) = x(t + T) — x(t). 

Fig. 4: Analogous specific heat C q of the linear (upper curve) and non-linear (lower 
curve) LD derived using Eq.(^). 
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